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We derive a new exact self-consistent crystalline condensate in the 1 + 1 dimensional chiral Gross- 
Neveu model. This also yields a new exact crystalline solution for the one dimensional Bogoliubov-de 
Gennes equations and the Eilenberger equation of semiclassical superconductivity. We show that the 
functional gap equation can be reduced to a solvable nonlinear equation, and discuss implications 
for the temperature-chemical potential phase diagram. 
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Interacting fermion systems describe a wide range of 
physical phenomena, from particle physics, to solid state 
and atomic physics [H, 0, M 0] ■ Important paradigms 
include the Peierls-Frohlich model of conduction [1, the 
Gorkov-Bogoliubov-de Gennes approach to superconduc- 
tivity [6j, and the Nambu-Jona Lasinio (NJL) model 
of symmetry breaking in particle physics pj. A 1 + 1- 
dimensional version of the NJL model, the NJL2 model 
[also known as the chiral Gross-Neveu model, XGN2] has 
been widely studied as it exhibits asymptotic freedom, 
dynamical mass generation, and chiral symmetry break- 
ing [H, [tj [l(| [lH . Surprisingly, the temperature-density 
phase diagram of this system is not yet fully understood. 
A gap equation analysis based on a homogeneous con- 
densate suggests its phase diagram is the same as its 
discrete-chiral cousin, the original Gross-Neveu (GN2) 
model |8|], while recent work finds an inhomogeneous 
Larkin-Ovchinikov-Fulde-Ferrell (LOFF) helical complex 
condensate ( "chiral spiral" ) below a critical temperature 
[T2I ]. In fact, the phase diagram of the discrete-chiral ver- 
sion, the GN2 model, has only recently been solved in the 
particle physics literature, analytically and on the 
lattice [14(. There is a crystalline phase at low temper- 
ature and high density, and this phase is characterized 
by a periodically inhomogeneous (real) condensate that 
solves exactly the gap equation. This phase is not seen 
in the old phase diagram which was based on a uniform 
condensate [llj]. Interestingly, this discrete-chiral GN2 
model (with vanishing bare fermion mass) is mathemat- 
ically equivalent to several models in condensed matter 
physics: the real periodic condensate may be identified 
with a polaron crystal in conducting polymers 0, [Til . 17 1 , 
with a p eriodic pair potential in quasi ID superconduc- 
tors 18, 13], and with the real order parameter for super- 
conductors in a ferromagnetic field [201 ] . This system also 
is a paradigm of the phenomenon of fermion number frac- 
tionalization [21, 22 1. Variants of such models also apply 
to ultracold fermionic systems, for which there are in- 
teresting new theoretical and experimental developments 

In this paper we present an analogous complex crys- 
talline condensate for the chiral GN system, the NJL2 
model. This condensate is an exact inhomogeneous so- 
lution to the gap equation, and also provides a new self- 



consistent solution to the Bogoliubov-de Gennes (BdG) 
Q and Eilenberger [25| equations of superconductivity. 
Our solution may also be relevant for chiral supercon- 
ductors and for incommensurate charge density waves in 
quasi ID systems [261 ] . which have chiral symmetry and 
an inherently complex order parameter. 

Consider the massless NJL2 model with Lagrangian 
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We have suppressed summation over N flavors, which 
makes the semiclassical gap equation analysis exact in the 
N — > 00 limit, a limit in which we can consistently discuss 
chiral symmetry breaking in 2D. The original GN2 model 
8] , without the pseudoscalar interaction term {xpi^ip^ , 
has a discrete chiral symmetry ip —> 7 5 ^- There are two 
equivalent ways to find self-consistent static condensates. 
First, introduce bosonic condensate fields, S = ipip and 
P = ■pi^tp, which we combine into a complex condensate 
field: A = S — i P = Me 1 ". Integrating out the fermion 
fields we obtain an effective action for the condensate as 
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The corresponding (complex) gap equation is 
6 



(2) 



A(x) = -2iNg" 



SA(x)* 
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i p - M{x) e~ 11 x{x) 
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One of the main results of this paper is that this gap 
equation can be reduced in an elementary manner to an 
ordinary differential equation, which moreover is solu- 
ble. A second approach to finding a self-consistent con- 
densate is to solve the relativistic Hartree-Fock prob- 
lem Hip = Eip, subject to the consistency condition 
{pip} — i(ipij 5 ip) — —A/g 2 , with single-particle Hamil- 
tonian 



A*(z) 



" dx 



H = -ij 5 — + 7 M(x) e~ 47 x{x) 
dx 

(4) 

This is the Bogoliubov-de Gennes (BdG), or Andreev, 
hamiltonian, with A (a;) playing the role of the order pa- 
rameter We have chosen Dirac matrices 7 = <ri, 
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7 1 = — i&2, and 7 5 = 03, to emphasize the natural com- 
plex combination A = S — iP. The key object in our 
analysis is the Gork'ov Green's function, and in particu- 
lar its coincident-point limit, the "diagonal resolvent": 



R(x;E) 
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This is a 2 x 2 matrix, and the spectral function is 
p(E) = ilmTro^ R(x; E + ie), where the trace is both 
a Dirac and spatial trace. Approximation methods, such 
as the gradient and semiclassical expansions, of R(x; E) 
have been widely studied [27j |. In one spatial dimension, 
R(x; E) can be written in terms of two independent solu- 
tions "01,2 to the Dirac/BdG equation Hip = Eip. That 
is, R(x; E) = (tpxtpj + V^^f) ci/(2W), with Wronskian 
W = —iip'[(J2'ip2- It follows immediately that R(x;E) 
satisfies the following first order equation: 

i?' 7 5 =i[ 7 5 (^- 7 °Me- l75x ),i?7 5 ] (6) 

In superconductivity, ([6]) is known as the Eilenberger 
equation 12511 . and in mathematical physics as the Dik'ii 
equation [28| |. It is also straightforward to show that 
R = R\ deti? = 1/4, and tr(i? 7 5 ) = 0. 

Our main observation is that the gap equation mo- 
tivates a self-consistent ansatz form of the 2x2 matrix 
R(x; E), and when this is combined with the identity ©, 
the exact self-consistent condensate and associated resol- 
vent are completely determined. To see this, note that 
the gap equation (J3]) can be viewed in two ways. First, 
write the log det term as 
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All dependence on A(x) resides in the spectral function 
p(E), via Tri?, and so the simplest solution to the gap 
equation (|3|) is for diagonal entries of the 2x2 matrix 
R [recall they are equal since tr(i? 7 5 ) = 0] to be linear 
in |A| 2 . On the other hand, we can also express the gap 
equation, by performing the functional derivatives, as 



A(x) = iNg 2 tr DiE [( 7 ° + 7 1 ) R(x; E)] 



(8) 



This suggests that the off-diagonal entries of R be pro- 
portional to A and A*. In fact, consistency between ([6]) 
and ((8|) introduces derivative terms, leading to 



R(x; E) = N{E) 



a(E) 
b(E)A' 
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The gap equation is satisfied since we find %teM{E) = 0. 
With this ansatz, the diagonal entries of ([6]) are identi- 
cally satisfied, while the off-diagonal entries imply that 
the condensate A satisfy the equation 

A" - 2|A| 2 A + % (b - 2E) A' - 2 (a - Eb) A = (10) 



This nonlinear Schrodinger equation (NLSE) is analyti- 
cally soluble, and all previously known examples of self- 
consistent condensates in GN2 and NJL2 are special 
cases. Furthermore, the corresponding Dirac/BdG equa- 
tion Hip — Eip is also exactly soluble, and has a spectrum 
consisting of a single band in the gap [13] • This provides 
an elementary explanation of the result from inverse scat- 
tering [iJElJEI that self-consistent ground state conden- 
sates are reflectionless (or finite-gap) with a single state 
(or band). For example, if we specialize to a real conden- 
sate, as is relevant for the [discrete-chiral] GN2 model, the 
[rescaled] self-consistent solution is A(x) = \fv sn(x; v) , 
which satisfies A" — 2 A 3 = — (1 + v)A : where sn is a Ja- 
cobi elliptic function, with elliptic parameter < v < 1 
13, 12 , HI, E3, H3] ■ The fermion spectrum has a single 
band in the gap, centered on E — 0, reflecting the charge- 
conjugation symmetry of the GN2 system. As v — > 1, the 
period becomes infinite and we obtain the famous kink, 
A(a;) = tanh(x), with a single bound state at E = 0. 
This mid-gap zero mode has many interesting physical 
consequences in polymer systems, and is the paradigm 
of the fractional fermion number phenomenon [2l( . It 
is also worth noting that for the [discrete-chiral] GN2 
model with a bare fermion mass, the NLSE (fT0|) acquires 
an inhomogeneous term, in which case the general solu- 
tion is written as A(x) = (("/) + C( x ) ~ C( x + 7)7 where 
C is the Weierstrass zeta function, and this represents a 
kink-antikink crystal (or bipolaron crystal in the polymer 
language [29(), which is a periodic generalization [l3( of 
the Dashen-Hasslacher-Neveu (DHN) kink-antikink solu- 
tion [9j: A(x) = coth(6) + tanh(x) — tanh(x + b). 

All these are real condensates, and are well-known. 
The only previously known complex condensates are (i) 
the simple "chiral spiral" or LOFF solution A(x) = 
Ae iqx ] and (ii) Shei's "twisted kink" solution [13], which 
we can express in complex form 



A{x) 
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9 is the angle through which the phase of the condensate 
rotates as x goes from —00 to +00. The single-particle 
fermion spectrum has a single bound state within the 
gap, located at E — cos (|) , as shown by the dashed line 
in Figure [TJ Consistency with the gap equation requires 
vanishing of the coefficient of A', which places condi- 
tions on parameters of the solution. For example, for 
Shei's solution the condition is that 0/(2n) is equal to 
the filling fraction of the bound level 13, 11 1 ■ We point 
out a simple physical interpretation of this condition in 
terms of conserved currents (see also 3l|). Both NJL2 
and GN2 models have a conserved current = ipj^ip. 
Since (j^ix)) — —iNtr DtE [l°"/ tJ 'R{x)], we see that, for 
a static condensate, current conservation follows trivially 
since (j 1 (x)) = 0; this is the physical origin of the identity 
trD (j 5 R) = 0. The NJL2 model also has a conserved ax- 
ial current, j£ = ipj^j^ip, and the Eilenberger equation 
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FIG. 1: Band spectrum of the BdG Hamiltonian Q for the 
complex crystal condensate (|13[) . as a function of the twist 
parameter 0. The dashed line denotes the infinite period limit 
[in this case v — > 0], which is Shei's solution (JTTJ - At 6 = 2K 
we recover the symmetric spectrum of the real kink crystal, 
relevant to the GN2 system. 



([6]) implies 

8^) = -iNtr(R') = 2S{x)$ifi>) - 2P(x)(^) (12) 

Axial current conservation then follows from the gap 
equation, using precisely the same condition on the coef- 
ficient of A'. 

We have found a nontrivial crystalline solution to (fT0|) 



A 



a{iAx + 8/2-K) 



(13) 



a(iAx-K)a{9/2) 
x exp [iA x (-C(0/2) + cs(0/2)) + 6*77/2 + i am(0/2)] 

Here, a and £ are Weierstrass sigma and zeta functions, 
A(6) = 2sd(0/4)cn(0/4), and r/ = ((K), with K(u) 
the elliptic half-period. For this condensate, both the 
amplitude and the phase are rr-dependent, as shown in 
Figs. [H [21 The essential parameters of the solution 
(fT3|) are: (i) the parameter 8 which [via (fl~5j) below] 
characterizes the chiral twist of A(x) over one period; 
and (ii) the elliptic parameter v which, together with 9, 
determines the crystal period. The spinor solutions of 
the Dirac/BdG equation can also be expressed explicitly 
in terms of elliptic functions, and one can perform the 
Hartree-Fock analysis, as in the GN2 system, to prove 
self-consistency of the crystalline condensate 30]. The 
spectrum is that of a Dirac particle with a single band 
in the gap, as shown in Fig. [TJ However, unlike the 
real case where the band lies symmetrically in the cen- 
ter of the gap, here the band is offset. Indeed, the band 
edges are given by: E\ = —1, E 2 = — 1 + 2 cn 2 (#/4; v), 
E 3 = -1 + 2cd 2 (6»/4; u), E 4 = +1, as shown in Fig. Q] 
This spectrum is a band version of the Shei spectrum, 
reducing to the Shei solution in the infinite period limit. 
It is also a deformation of the kink crystal spectrum of 
the [discrete-chiral] GN2 system, reducing to that case 
when 9 = 2K. We find the exact diagonal resolvent ((9]) 
with 
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1 



*y/{E*-l)(E-E 2 )(E-E 3 ) 



(14) 



b(E) = 2E - (E 2 + E 3 ), and a(E) = 2(E + l)(E - 1 - 
(E 2 + E 3 )/2) + l + (E 2 + E 3 ) - (E 2 - E 3 ) 2 /A. 

Under a shift through one period L = 2K.'/A of the 
crystal, the BdG Hamiltonian is invariant up to a global 
chiral rotation through an angle (p: 



H{x 



L) 



e n v H{x)t 
K'[-C(0/2) 



-cs(6»/2) 



(15) 

ieC(iK')/(2K')] 



The solutions to Hip = Eip acquire a chiral rotation 
and a Bloch phase under a period shift, ip(x + L) = 
e ikL e itp-y jp( x y xhe [relativistic] Bloch momentum k is 
related to the the spectral function by p{E) = dk/dE. 



M(x) (Solid) 




X(x) (Dashed) 



FIG. 2: The amplitude M(x) [solid line] and phase x(x) 
[dashed line] of A = Me lx in (]13[l. over several periods. The 
amplitude is periodic while the phase rotates by an angle 2tp 
each period. 




FIG. 3: Crystalline complex kink fj 13p plotted as a func- 
tion of x over one period. The cross-section denotes A(x) = 
S(x) — iP(x), and indicates a net rotation through the twist 
parameter 2ip over one period. 

To conclude, we discuss briefly the implications of this 
self-consistent solution to the gap equation for the (T, fi) 
phase diagram of the NJL2 model. Recall that the cor- 
responding real condensate characterizes the inhomoge- 
neous crystalline phase of the GN2 system [13| . For the 
NJL 2 model, the Ginzburg-Landau (GL) approach shows 
that the " chiral spiral" phase identified in [12( has a richer 
structure, characterized by a crystalline complex conden- 
sate of the form (fTB")) . Near the tricritical point of the 
massless NJL2 system the GL effective Lagrangian is 

C GL = C o + C 2|A| 2 + c 3 Im[A(A')1+c 4 [|A| 4 + |A 
+c 5 Im[(A"-3|A| 2 A) (A')*] 
+c 6 [2|A| 6 + 8|A| 2 |A'| 2 + 2Re ((A') 2 (A*) 2 ; 
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(16) 
A"| 2 ] 
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Here the coefficients c„ are known functions of T and fi 
[32 , l33fl . (This GL approach has been used previously 
in 33j to describe the phase diagram of the massive and 
massless NJL 2 models, in the vicinity of the tricritical 
point. For the massive NJL2 model no (complex) ex- 
act solution to the gap equation is known, so [33] is the 
current state-of-the-art for the massive model.) In the 
GN2 model, which has a real condensate, there is a tri- 
critical point at T = 0.3183, fi — 0.6082; given by the 
point C2 = C4 = 0. In the [chiral] NJL2 model, which 
has a complex condensate, there is a tricritical point at 
T = 0.5669, /1 = 0; given by C2 = C3 = 0. To search for 
possible crystalline phases near the tricritical point we 
keep terms up to C4, and study the effective equation of 
motion 



c 4 A" - ic 3 A' - (c 2 + 2c 4 |A| 2 ) A = 



(17) 



Note that this has precisely the same form as the NLSE 
(j 1 0[) found from the gap equation and the Eilenberger 
equation. Thus, we can use our solution (|13p as a 
variational ansatz and compute the free energy. We 
have found that the resulting free energy is lower than 
that of the LOFF-form "chiral spiral" variational ansatz 
A(x) = Ae iqx , for T < T c . Noting that our ansatz re- 
duces to the chiral spiral in the perturbative limit (just 
as the GN 2 self-consistent crystal condensate reduces to 
the LOFF form A(x) = Asin(qx) in this limit [3]), we 
are led to the phase diagram shown in Fig. 0] For a de- 
tailed discussion of the phase diagram, based on the full 
free energy (beyond the GL approximation), using our 
exact spectral data, see [301 ] . 



T 



LOFF "chiral spiral" 
twisted kink crystal 



FIG. 4: Phase diagram of the chiral NJL2 model, from a 
Ginzburg-Landau analysis based on a crystalline condensate 
of the form in (|13[1 . Along the T axis (below T c ) the conden- 
sate is of the form of Shei's twisted kink while on the 
T = T c line the condensate reduces to the form of the LOFF 
chiral spiral The condensate is depicted as in Fig. [3] 
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